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Abstract. This work lies across three areas (in the title) of investigation that are by themselves of 
independent interest. A problem that arose in quantum computing led us to a link that tied these 
areas together. This link consists of a single formal power series with a multifaced interpretation. 
The deeper exploration of this link yielded results as well as methods for solving some numerical 
problems in each of these separate areas. 
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1. Introduction 



Since our work may be of interest to audiences of varied background we will try to keep our notation 
as elementary as possible and entirely self contained. 

The problem in invariant theory that was the point of departure in our investigation is best stated 



version. Given two matrices A = 


"an 


012' 


and B = 


bii 612" 




-021 


022- 




-^21 ^22- 



of determinants 1, or equivalently in S'L[2] := 5L(2,C), we recall that their tensor product may be 
written in the block form 



A<S)B 



aiiB ai2B 
021-6 022-6 



(1) 



We also recall that the action of a matrix M — [mijjfj^i on a polynomial P{x) in R„ 
C[a;i, a;2, . . . , x„] may be defined by setting 



TmP{x) = P{xM), 



(2) 



where the symbol xM is to be interpreted as multiplication of a row n-vector by an n x n matrix. This 
given, we denote by j^^^[2]'8'S^[2] Tcing of polynomials in R4 that are invariant under the action of 
A (g) -B for all pairs A, B e SL[2]. In symbols 



R 



S-L[2]«iS-L[2] 



{PeR4 : Ta^bP{x) ^ P{x)}. 



(3) 



Since the action in ^ preserves degree and homogeneity, R4 
it decomposes into the direct sum 



SL12]»SL12] . 



is graded, and as a vector space 



R 



SL[2]0SL[2] _ ^^Jf ^j^S-L[2]®SL[2]^ 

m>0 



where the m^^ direct summand here denotes the subspace consisting of the SL[2] (g) 5-L[2]-invariants 
that are homogenous of degree to. The natural problem then arises to determine the Hilbcrt series 



W2iq) = g"dim7^„(R; 



SL[2](»SL12] 



m>0 
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Now note that usmg ([T]) iteratively we can define the fc-fold tensor product Ai ^ A2 (E) ■ ■ ■ (E) Ak, and 
thus extend ^ to its general form 

and set 

m>Q 

Remarkably, to this date only the series W2{q),W3{q),W4{q),W5{q) are known explicitly. Moreover, 
although the three series W2{q),W3{q),W4{q) may be hand computed, so far W5{q) has only been 
obtained by computer. 

The third named author, using branching tables calculated in [8], was able to predict the explicit 
form of W^^q) by computing a sufficient number of its coefficients. The computation of these tables 
took approximately 50 hours using an array of 9 computers. 

The series W4{q), ^5(5) first appeared in print in works of Luque-Thibon [5], [B] which were moti- 
vated by the same problem of quantum computing. We understand that their computation of ^5(5) 
was carried out by a brute force use of the partial fraction algorithm of the fourth named author, and 
it required several hours with the computers of that time. 

The present work was carried out whilst unaware of the work of Luque-Thibon. Our main goal is 
to acquire a theoretical understanding of the combinatorics underlying such Hilbert series and give a 
more direct construction of W5{q) and perhaps bring We{q) within reach of present computers. 

Fortunately, as is often the case with a difficult problem, the methods that are developed to solve 
it may be more significant than the problem itself. This is no exception as we shall see. 

Let us recall that the pointwise product of two characters x*-^-* and x^^^ of the symmetric group Sn is 
also a character of Sn, and we shall denote it here by x^^-* ©x''^'- This is usually called the Kronecker 
product of x^^'' and x^^^- An outstanding yet unsolved problem is to obtain a combinatorial rule for 
the computation of the integer 

'^A(i>,A(2),...,A('=) i^) 

giving the multiplicity of x^ in the Kronecker product x^* ' © X"*** '©■••© X^* ' ■ Here x"^ and each 
X"^* ' are irreducible Young characters of Sn- Using the Frobenius map F that sends the irreducible 
character x^ onto the Schur function S\, we can define the Kronecker product of two homogeneous 
symmetric functions of the same degree / and g by setting 

/0g = F((F-i/)0(F-\9)). 

With this notation the coefficient in ^ may also be written in the form 

'^A<i),A<2),...,A<'=) = (*A(i) © Sx(2) © ■ ■ • © Sx(k) , Sa), 

where ( , ) denotes the customary Hall scalar product of symmetric polynomials. The relevancy of 
all this to the previous problem is a consequence of the following identity. 

Theorem 1.1. 

Wkiq) = q^'^{sd,d © Sd,d © • • • © SdA , S2d) (5) 

where, in each term, the Kronecker product has k factors. 

For this reason, we will often refer to the task of constructing Wk{q) as the Sdd Problem. Using this 
connection and some auxiliary results on the Kronecker product of symmetric functions we derived in 
[3] that 

W2{q) = -^, Wsiq)^-^, W4{q) = - ^ ^. (6) 

1-q^ [1- q^){l- q-^y{l- q^) 
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Although this approach is worth pursuing (see [3]), the present investigation led us to another sur- 
prising facet of this problem. 

Let us start with a special case. We are asked to place integer weights on the vertices of the unit 
square so that all the sides have equal weights. Denoting by Pqo, Pqi, Pio, Pii the vertices (see figure) 
and by pqq, poij PiOi Pii their corresponding weights, we are led to the following Diophantinc system. 



Poo+Pai - Pio-Pii = 
Poo -Pai + Pio -Pii = 0' 



The general solution to this problem may be expressed as the formal series 

1 




F2iyoo,yoi,yio,yii) 



E POO POl PlO Pll 
yoo yoi tfio yii 

peS2 



(1 - yoo2/ii)(l - yoi2/io) 



In particular, making the substitution yoo = j/oi = yio ~ yii = q we derive that the enumerator of 
solutions by total weight is given by the generating function 

G2(g) = ^771,(2)92" 1 

d>0 



(1 



with md{2) giving the number of solutions of total weight 2d. 

This problem generalizes to arbitrary dimensions. That is we seek to enumerate the distinct ways 
of placing weights on the vertices of the unit fc-dimensional hypercube so that all hyperfaces have the 
same weight. Denoting by pe^e^.-.e^ the weight we place on the vertex of coordinates (ei, £2, • ■ • , Cfc) we 
obtain a Diophantine system Sk of k equations in the 2*^ variables {?5eie2 - efc}ei=o,i- 
For instance, using this notation, for the 3-dimensional cube we obtain the system 





Pooo + Pool + POlO + Poll - PlOO - Pioi - PllO - Pill 



Pooo + Pool - POlO - Poll + PlOO + PlOl 
POOO - Pool + POlO - Poll + PlOO - PlOl 

In this case the enumerator of solutions by total weight is 

1 



PllO 
PllO 



Pill 
Pill 



d>0 



md{i)q 



2d 



{l-q^ni-q^f 



The relevance of all this to the previous problem is a consequence of the following identity. 

Theorem 1.2. Denoting by md{k) the number of solutions of the system Sk of total weight 2d and 
setting 

(7) 



d>0 



Gk{q) = q^'^{hd,d © hd,d • • • hd,d , 5*2^), 



d>0 

where, hd^d denotes the homogenous basis element indexed by the two part partition {d,d), and in each 
term, the Kronecker product has k factors. 

For this reason, we will refer to the task of constructing the series Gk{q) as the Hdd Problem. 

Theorem 11.21 shows that the algorithmic machinery of Diophantine analysis may be used in the 
construction of generating functions of Kronecker coefficients as well as Hilbert series of ring of in- 
variants. More precisely we are referring here to the constant term methods of MacMahon partition 
analysis which have been recently translated into computer software by Andrews et al. [1] and Xin 

m. 
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To see what this leads to, we start by noting that using MacMahon's approach the solutions of ^2 
may be obtained by the following identity 



F2(yoo,2/oi,2/io,2/ii)= E E E J2yooy'ory'ioy'il'^'i 



Poo+Pol — Pio— Pii _Poo— POl +P10— Pll 



as 



Poo>0poi>0pio>0pii>0 

where the symbol " o o" denotes the operator of taking the constant term in 01,02- This identity 
may also be written in the form 

1 



F2{yoo,yoi,yw,yii) 



(1 - yooaia2)(l - 2/oiai/a2)(l - 2;ioa2/ai)(l - 2/11/0102) 
In particular the enumerator of the solutions of ^2 by total weight may be computed from the identity 

G2{q) ^ 



More generally we have 



(1 - gai02)(l - 9ai/o2)(l - go2/oi)(l - g/0102) 



n (i-9n.,esa'/nj^s«j 

SC[l,fc] 



(8) 



a\'a);---a". 



where we use (and will often use) [m, n] to denote the set { to, to + 1, . . . , n }. Now, standard methods 
of Invariant Theory yield that we also have 



Wk{q) 



nti(i-< 



SC[l,fe] 



(9) 



A comparison of ^ and ^ strongly suggests that a close study of the combinatorics of Diophantine 
systems such as Sk should yield a more revealing path to the construction of such Hilbert series. This 
idea turned out to be fruitful, as we shall see, in that it permitted the solution of a variety of similar 
problems (see [3], [1]). In particular, we were eventually able to obtain that 



(1 - 9)9(1 - g2)8(l _ g3)6(l _ g4)3(i _ ^5) ' 



(10) 



with 



N5 = q^^ + 7q^^ + 220q*2 + 26069"*^ + 242299"*" + 169840q^9 + 951944q3* 



4391259g3^ + 17128360g'''' + 57582491?^'' + 169556652g^^ + 442817680? 



^36 



35 



^34 



^33 



+ 1036416952g32 + 21921916075^* + 4219669696?^° + 7433573145g29 + 120413052719^8 
+ 18003453305g^^ + 2492175141692^^ + 32017113319?^^ + 382432748519^'' + 425248150139^^ 
+ 440524404329^2 _^ 425248150139^* + 3824327485192" + 320171133199*^ + 249217514169*^ 
+ 180034533059*^ + 120413052719*^ + 74335731459*^ + 42196696969*"* + 21921916079*^ 
+ 10364169529*2 + 4428176809** + 1695566529*" + 575824919^ + 171283609^ + 43912599^ 
+ 9519449^ + 1698409^ + 242299"* + 26069^ + 22O92 + 79 + 1. 

Surprisingly, the presence of the numerator factor in ([9]) absent in ([8]) does not increase the complexity 
of the result, as we see by comparing (jlO|) with the Luque-Thibon result 



W^5(V9) 



{l-q^Y{l-q^){l^q^f{l-q^)[l-q^r 
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with 

Pg = + + 16g^° + 9/^ + 98g4^ + 154g4^ + 465g^^ + 915^''^ + 2042g4* + 37949*3 _^ Y263q'^^ 

+ 126889**^ + 21198g*° + 34323^^^ + 52205?^^ + 77068g3^ + 108458g3^ + 1474239^^ + 191794g34 
+ 241863g33 + 2926899^2 + 3A2207q^^ + 386980g3° + 421057?^^ + 443990g28 + 451398g2^ 
+ 443990?^^ + 421057g^^ + 3869809^"* + 3422075^^ + 292689g^^ + 241863g2^ + 191794(7^° 
+ 147423g^^ + 108458g^^ + 77068g" + 52205g^^ + 34323g^^ + 21198g" + 12688g" 
+ 7263^12 + 3794^" + 2042gi° + 9159^ + 465g^ + 1549^ + 98g^ + 9q^ + 16q^ + q^ + l. 

It should be apparent from the size of the numerators of W^iq) and G5{q) that the problem of 
computing these rational functions explodes beyond k — i. In fact it develops that all available 
computer packages (including Omega and Latte) fail to directly compute the constant terms in ([8]) for 
k = 5. This notwithstanding, we were eventually able to get the partial fraction algorithm of Xin [10| 
to deliver us G5(g). 

This paper covers the variety of techniques we developed in our efforts to compute these remarkable 
rational functions. Our efforts in obtaining We{q) and GQ{q) are still in progress, so far they only 
resulted in reducing the computer time required to obtain W5(<?) and G5{q). Using combinatorial 
ideas, group actions, in conjunction with the partial fraction algorithm of Xin, we developed three 
essentially distinct algorithms for computing these rational functions as well as other closely related 
families. Our most successful algorithm reduces the computation time for W^^q) down to about five 
minutes. The crucial feature of this algorithm is an inductive process for successively computing the 
series Gk{q) and Wk{q), based on a surprising role of divided differences. 

This paper is the extended version of [5]. We organize the contents in 5 sections. Section 1 is this 
introduction. In Section [5] we relate these Hilbert series to constant terms and derive a collection of 
identities to be used in later sections. In Section [3] we develop the combinatorial model that reduces 
the computation of our Kronecker products to solutions of Diophantinc systems. In Section [4] we 
develop the divided difference algorithm for the computation of the complete generating functions 
yielding Wk{q) and Gk{q)- In Section [U after an illustration of what can be done with bare hands 
we expand the combinatorial ideas acquired from this experimentation into our three algorithms that 
yielded G^lq) and our fastest computation of ^5(9). 

The readers are referred to the papers of Luque-Thibon [5] , [5] and Wallach [H] for an understanding 
of how these Hilbert series are related to problem arising in the study of quantum computing. 



2. Hilbert series of invariants as constant terms 

Let us recall that given two matrices A = [fly]™^]^ and B = [bij]2j^i we use the notation B 
to denote the nm x nm block matrix A^ B = [a^ =1 • For instance, if m = n = 2, then 

012^*12 

ai2&22 

022^12 

a22&22_ 

Here and in the following, we define TaP{x) to be the action of an m x m matrix A ~ [oij]™ =1 on a 
polynomial P{x) = P{xi,X2, ■ ■ ■ ,Xm) hi := C[xi,X2, ■ ■ ■ by 

m m m 

TaP{xi,X2, ■ ■ . jim) = p(^^^Xiaii , ^Xiai2 , . ■ . , (11) 

In matrix notation (viewing x = (xi, a;2, . . . , Xm) as a row vector) we may simply rewrite this as 

TAPix)=P{xA). 



A(E)B = 



anbii 


aii6i2 


012611 


aii&2i 


011622 


012621 


021611 


021612 


022611 


021621 


021622 


O22621 
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Recall that if G is a group of to x to matrices we say that P is G-invariant if and only if 

TaP{x) = P{x) V A EG. 

The subspace of of G-invariant polynomials is usually denoted R^. Clearly, the action in pT|) 
preserves homogeneity and degree. Thus we have the direct sum decomposition 

R^ = 7^o(R^0 © Wi(R^) © W2(R^) © • • • © HdiRg) © • • • 

where Tid^Rm) denotes the subspace of G-invariants that arc homogeneous of degree d. The Hilbert 
series of R^ is simply given by the formal power series 

FG(g) -^g''dim(Hd(R^; 

This is a well defined formal power series since dimTid^Rm) — "^ii^^ ^^d(R-n 

When G is a finite group the Hilbert series Fciq) is immediately obtained from Molien's formula 



/d+m-l\ 
\ m-1 / 



For an infinite group G which possess a unit invariant measure cu this identity becomes 

PG{q)= f j-TT^ T^du;. (12) 

JaeG det (/ - qA) 

For the present developments we need to specialize all this to the case G = S'L[2]®'^, that is the group 
of 2*^ X 2*^ matrices obtained by tensoring a fc-tuple of elements of 5i[2]. More precisely 

SL[2f'^'' ^ {Ai<»A2<»---®Ak : A,e SL[2] V i ^ 1,2, . . . ,k}. (13) 

Our first task in this section is to derive the identity in That is 

Theorem 2.1. Setting for k > 1 

Wk{q) = FsLl2]^.{q) = ^g-^dim (HdCEl'^t^'^'"')) , (14) 
d>0 



SC[l,fe] 



(15) 



avdo ■ ■ - a. 



We need the following result. 



Proposition 2.2. IfQ(ai, 02, . . . , a^) is a Laurent polynomial in C[ai, 02, . . . , ttk', 1/ai, l/a2, ■ • ■ , l/ofe] 
then 

^) r ••■ r Q(e''Se*'%...,e^''=)d0id02---d0fe = Q(ai,a2,...,afc) (16) 

Proof. By multilinearity. it suffices to consider Q{ai, 02, . . . , ak) = 0,1^^0,2' ■ ■ ■ a^^ , in which case psp 
obviously holds. □ 



Proof of Theorem \2.1\ To keep our exposition within reasonable limits we will need to assume here 
some well known facts (see [5] for proofs). Since SL[2\ has no finite measure the first step is to note 
that a polynomial P{x) E C[xi,X2, ■ ■ ■ , X2k] is S'L [2]® '^-invariant if and only if it is S'f7[2]®''-invariant, 
where SU[2] := S'[/(2, C) and as in ([E 



SU[2f'' ^ {Ai®A2®---^Ak : A,E SU[2] V i = 1, 2, . . . , fc}. 
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In particular we derive that -Fsl[2]»'= (?) = ^sc/[2]»'' (?)• This fact aUows us to compute ^5L[2]»fc (9) 
using MoUen's identity (|12p. Note however that if 

Ai(E)A2(E)---(E)Ak 

and Ai has eigenvahies i,; , 1 /U then (using plethistic notation) we have 

Denoting by dcjj the invariant measure of the i^^ copy of S'[/[2] we see that (fT2|) reduces to 

FsUl2]e>-{<l) = l"l / /lm[(tl + lAl)---(tA; + lAfe)]dt-'l---rf'^fe. 

Now it is well know that if an integrand f{A) of S'L/[2] is invariant under conjugation then 



/ f{A)dw ^ - r f 

JSU\2] ^ J -TV 



e^" 




sin^ Ode. 



(17) 



This identity converts the right-hand side of p7|) to 

^ • • + • • ■ {^"'' + e"*''")] sin^ O^--- sin^ 0,, • • • (18) 



The substitution 

reduces the coefficient of to 
1 



sin^ 6*, 



J — TT ^ — TT 



) dei---d0k. (19) 



However the factor /im[(e'^^ + e '^1 ) • • • (6'^*= + e '^'=)] is invariant under any of the interchanges 
e'^j < — > e^*^J . Thus the integral in (fTOl) may be simplified to 



(2^) 



l^f---fh^ [{e^'^ + e-^'^ )■■■ {e^'^ + e-«^ )] ^ (l - e^^'^ )de,--- d0. 



Proposition [2T2] then yields that this integral may be computed as the constant term 

k 

hm [{ai + l/ai)(a2 + l/aa) • • • (a^ + l/ofc)] H ~ ' 



Using this in we derive that 

Fsumax- (q) = XI '?"'*™ [("1 + l/ai)(a2 + 1/02) • • • (ofe + l/ofc)] H " 

m>0 



i>0 



n 



This completes the proof of Theorem 12.11 



i=l 

n(i- 



□ 
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Note that if we restrict our action of SU[2]^^ to the subgroup of matrices 
then a similar use of Mohen's theorem yields the following result. 

Theorem 2.3. The Hilbert series of the ring of invariants R^i? is given by the constant term 

1 



= (<?) = 



SC[l,fc] 



(20) 



Proof. The integrand 1/ det(l — qA) is the same as in the previous proof and only the Haar measure 
changes. In this case we must take dw ~ d9id92 ■ ■ ■ dOk/i^-K)^ in and Molien's theorem gives 



1 



Thus follows from Proposition 



SC[l,fc] 



deide2 ■ ■ ■ dOk 



□ 



Remark 2.4. There is another path leading to the same result that is worth mentioning here since 
it gives a direct way of connecting Invariants to Diophantine systems. For notational simplicity we 
will deal with the case k = 2>. Note that the clement 



"ti 0' 




h 


0' 




'ts 0' 


ti 


® 





^2. 




^3 



is none other than the 8x8 diagonal matrix 



A(ti,t2,t3) 



tlt2t'i 













tlt2/t3 













tlt3lt2 













iiA2i3 












t2t3/tl 

























tz/tit2 













lAii2i3 



This gives that for any monomial ^ ■ ■ ■ x^^ we have 

A{tl t2 ig) j;P = fPl+P2+P3+P4-P5-P6-p--PS ^Pl+P2-P3-P4+P5+P6-P7-P8 ^|l-P2+P3-P4+P5-P6+P7-PS ^^P 

Thus all the monomials arc eigenvectors and a polynomial P(xi,X2, . . . ,xs) will be invariant if and 
only if all its monomials are eigenvectors of eigenvalue 1. It then follows that the Hilbert scries Fr^^x (q) 



of C[xi,X2, 



is obtained by q-counting these monomials by total degree. That is g-counting 



by the statistic pi + P2 + Ps + P4 + P5 + Pe + Pr 

S3 

and MacMahon partition analysis gives 

Frpsia {q) 



the solutions of the Diophantine system 



Pi+P2+P3+Pi-P5-P6-P7-Ps = 

P1+P2-P3-P4+ Pb +P6-P7 -PS = 
Pi - P2+P3 - P4+P5 - P6 +P7 - PS = Q 



(21) 



l — qaia2a3 l — qa\a2/a^ l — qaia-3,/a2 l — qa\/a2a^ l — qa2a^/a\ 1 — 902/0103 l — qa^/a\a2 1 — qa/ aia^a^ 000' 

This gives another proof of the case A: = 3 of (PD)l . It is also clear that the same argument can be 
used for all > 3 as well. 
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Remark 2.5. Full information about the solutions of our systems is given by the complete generating 
function 



Fk{xi,X2, .. .,X2k) = ^ X]": 
peSk 



,P2 



P2k 



(22) 



Using the notation adopted for ^3 in (|2ip . our system Sk may be written in vector form 

PlVl + P2V2 + • • • + P2kV2k = 0, 

where Vi, V2, . . . , ^2'" ^re the fc-vectors (±1, ±1, . . . , ±1) yielding the vertices of the hypercube of 
semiside 1 centered at the origin. In this notation, MacMahon partition analysis gives that the 
rational function in (|22|) is obtained by taking the constant term 



Fk{xi,X2, 



1 X2'' ) 



with the Ai Laurent monomials in ai, 02, . . . afe which may be written in the form 

where £162 ■ ■ ■ are the binary digits of i — 1. 

In the same vein the companion rational function W{x\^ 2:2, ... , X2k) associated to the Sdd problem 
is obtained by taking the constant term 

2*.- 



Wk(x\,X2,. ■ ■ ,X2k) = 11(1 -"1)11^3" 

7 = 1 i=l ^ 



XiAi 



(23) 



a';a'i---a" 



Of course we have 

Gk{q) = Fk{xi,X2, ■ . • ,X2fc) 



and Wkiq) = Wk{xi,X2, ■ ■ ■ ,X2k) 



In Section [4] we will show that, at least in principle, these rational functions could be constructed by 
a succession of elementary steps interspersed by single constant term extractions. 

3. DiOPHANTINE SYSTEMS, CONSTANT TERMS AND KRONECKER PRODUCTS 

We have seen, by MacMahon partition analysis, that the generating function Gk{q) defined in ([7]), 
which counts solutions of the Diophantine system Sk, is given by the constant term identity in ([5]): 



Gkiq) 



SC[l,fc] 



(24) 



a"a"---a". 



In the last section we proved (in Theorem 12. ip that the Hilbert scries Wk{q) of invariants in (jl4p is 
given by the constant term 

nti(i-«?) 



Wk{q) 



SC[l,fc] 



(25) 



a';a^---a". 



A comparison of ((24)) and ((25)) clearly suggests that these two results must be connected. This 
connection has a beautiful combinatorial underpinning which leads to another interpretation of the 
these remarkable constant terms. The idea is best explained in the simplest case k = 2. Then 
reduces to 



W2iq) = 



(1 - gaia2)(l - qai/a2){l ~ ga2/ai)(l - q/aia2) 



10 
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Expanding the inner rational function as product of four formal power series in q we get 



^^w= E E E E'^' 

Poo>0 poi >0 pio>0 pii>0 



poo +Poi+Pio+Pi1qP"o+Poi -pill -pii^Poo-poi+pio-pii 



On 



E 

Poo>0 poi>0 pio>0 pii>0 



E E E 



+P01+P10+P11 „Poo+Poi-pio-pii+2 Poo-Poi+Pio-Pii 
"1 "2 



(26) 



E E E E 

Poo>0poi>0pio>0pii>0 



Poo+Poi+Pio+Pii ^Poo+Poi-pio-pii^Poo-poi+Pio- 
tt-i tin 



1+2 



^poo+poi+pio 



E 

Poo>0poi>0pio>0pii>0 



+P11 a' 



Poo+Poi-Pio-Pii+2^Poo-Poi+Pio-pii+2 



(27) 



Now by MacMahon partition analysis, the the i^^ term counts solutions of the Diophantine system 

Poo + Poi - Pw - Pii = Ci 
Poo - Poi + Pio - Pii = di 
where {ci,di) equals (0, 0), (— 2, 0), (0, — 2), (— 2, — 2) for i = 1,2,3,4, respectively. Note that the first 
term of ((26)) is none other than (|24p for fc = 2. 

Applying the same decomposition in the general case we see that the series Wk {q) may be viewed 
as the end product of an inclusion exclusion process applied to a family of Diophantine systems. 
To derive some further consequences of this fact, it is more convenient to use another combinatorial 
model for these systems. In this alternate model our family of objects consists of the collection JF^ of 
d-subsets of the 2(i-element set 

02d = {l,2,3,...,2d}. 

For a given A = {1 < ii < i2 < ■ ■ ■ < id l£ 2c?} e J-^ and a in the symmetric group 5*2^ we set 

(tA^ {ai^,ai2, . . . , ai^}. 
This clearly defines an action of S2d on J^d as well as on the fe-fold cartesian product 

Theorem 3.1. The number md{k) of solutions of the Diophantine system Sk is equal to the number 
of orbits in the action of S2d on . 

Proof. It will be sufficient to see this for k = 2. Then leaving d generic we can visualize an element 
of !Fd X !Fd by the Ven diagram of Figure [TJ There we have depicted the pair (^1,^2) as it lies in 
f22d- Using these two sets we can decompose il2d into 4 parts labeled by Aqo, ^01, ^10, ^ii- More 
precisely "Aqo" labels the set AiH A2, "Aqi" labels the set Ai n ''A2, "Aio" labels the set ""Ai n A2 
and "All" labels the set "^i n ''A2. Here we use '"!A,; " to denote the complement of Ai in ^l2d- This 
given, if we let PoOiPoiiPiOiPii denote the respective cardinalities of these sets, the condition that the 
pair (^1,^2) belongs to J-d x Td yields that we must have 

Poo + Poi + Pio + Pii = 2d 
Poo +P01 = l^il = d 
Poo +P10 = 1^2! = d 
Note that this system of equations is equivalent to the system 

Poo + Poi + Pio + Pii = 2d 
Poo +P01 -Pio -Pii = . 
Poo -Poi +P10 -pii = 

It is easily seen that for any solution (pooiPoiiPiOiPii) of this system, we can immediately construct 
a pair of subsets (^1,^2) & ^d x ^d by simply filling the sets ^00,^01)^10,^11 in the diagram of 
Figure [1] with pooiPoijPiOjPii respective elements from the set ^2^. Moreover, any two such fillings 
can be seen to be images of each other under suitable permutations of S2d- In other words by this 
construction we obtain a bijection between the orbits of J^d x J^d under S2d and the solutions of the 
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Figure 1. The Ven diagram for JTj. 



system 52 we have previously encountered. 
foUows by an entirely analogous argument. 



This proves the theorem for fc = 2. 



The general case 
□ 



Now we are ready to prove Theorem II .21 and then Theorcm ll.il 



Proof of Theorem \1.2l We are to show that 

mk{d) ^ {hd,d hd,d ■ ■ ■ hd.d , S2d)- 



(28) 



It is well known that a transitive action of a group G on a set f2 is equivalent to the action of G on 
the left G-cosets of the stabilizer of any element of ft. In our case, pick the subset [1, d] of 02;;. Then 



the stabilizer is the Young subgroup Si 



X Si 



[d+l,2d] 



of 52^ and thus the Frobenius characteristic 



of this action is the homogeneous basis element hd^d = hdhd- It follows then that the Frobenius 
characteristic of the action of 5*21; on the fc-tuples {Ai,A2, . . . ,Ak) of d-subsets of il2d is given by the 
fc-fold Kronecker product hd,d • • • hd,d- Therefore the scalar product 

{hd,d hd,d • • • hd,d , S2d) 

yields the multiplicity of the trivial under this action. But it is well known, and easy to see that this 
multiplicity is also equal to the number of orbits under this action. Thus (|28p follows by Theorem 

o □ 



Proof of Theorem ] 1 . 11 Again we will only need to do it for k ~ 2. To this end note that by Theorem 
11.21 the number of solutions of the system S2 in is given by the scalar product 

{hd,dQhd,d , S2d)- (29) 

In the same vein we see that the number of solutions to the system 5| in ((TT]) may be viewed as the 
number of orbits in the action of S2d on the pairs of subsets (74.1,^2) of n2d where |^2| = I "^^2! and 
I All = I'^Ail + 2. We have seen that the Frobenius characteristic of the action of S2d on subsets of 
cardinality d is hd^d- On the other hand the action of S2d on sets of cardinality d + 1 is equivalent to 
the action of S2d on left cosets of 5'[i,d+i] ^ S[d+2,2d] yielding that the Frobenius characteristic for this 
action is hd+ihd-i- Thus the Frobenius characteristic of the action of S2d on such pairs must be the 
Kronecker product 

hd+ihd-i hdhd- 

It then follows that the number of solutions of the system 5| is given by the scalar product 

{hd+ihd-i hdhd , S2d)- (30) 

The same reasoning gives that the number of solutions of the systems iS| and 5| in (j27p are given by 
the scalar products 



{hdhd hd+ihd-i , S2d) and {hd+ihd-i hd+ihd-i , S2d)- 
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It follows then that the coefficient of q^'^ in the alternating sum of formal power series in ((26|) is none 
other than the following alternating sum of the scalar products in (PU)) . ([50)1 and (pij) . 

W2{q) = (hdhd © hdhd , S2<i) - {hd+ihd~i /id/id , S2d) 

q2d 

- {hdhd hd+ihd-1 , S2d) + {hd+ihd-i hd+ihd-i , S2d) 
= {{hdhd - hd+ihd-i) (/id/id - hd+ihd-i) , S2d) = {sd.d Sd,d , S2d)- 
Summing over d gives 



d>0 



An entirely analogous argument proves the general identity in ([5]). 

4. Enter divided difference operators 



□ 



There is a truly remarkable approach to the solutions of a variety of constant term problems which 
exhibit the same types of symmetries of the Hdd and Sdd problems. We will introduce the approach 
in some simple cases first. We define the double of the Diophantine system 



52 



PI+P2- Pa - P4 = 
Pi - P2 + P3 - Pi = 



SS2 = 



to be the system 

Pi + P2 - P3 - Pi + P5 + P6 - P7 - P8 = 
P1-P2+ P3 -Pi + P5~P6+P7-P8^0' 

As we can easily see we have simply repeated twice each linear form and appropriately increased the 
indices of the variables. Now suppose that we are in possession of the complete generating function 
of that is 



X2 x'^ 



FS2{X1,X2,X3,X4) = X : 

P&S2 

We claim that the complete generating function of SS2 is simply given by 

^552(2:1, 2:2, ■ ■ ■ ,Xs) = Si^5S2fih,7(>i,sPs2ixi, X2, X3, X4^), 



(32) 



where for any pair of indices we let Sij denote the divided difference operator defined for any 
function f{x) by 

_ f{x)-f{x)\ 



Proof of ([32]) . By MacMahon partition analysis we have 

1 1 

^52 (2^1, 2:2, 2^3, 2:4) 



1 



(1 - 2:10102) (1 - a;2ai/a2) (1 - 0:302/01) (1 - 0:4/0102) 



Now note that since 
1 



(1 - xiai02) 
we obtain similarly 



1 



1 



1 



ai02 



(33) 



(1 - 2:10102) (1 — 2:50102) / 2:1 — 2:5 (1 — 2:iaia2)(l — 2:50102) 
1 ai/a2 



h.7 



(1 


- 2:201/02) 


(1- 


2:20i/02)(l - 


2:501/02) ' 




1 




02/01 




(1 


- 2:302/01) 


(1- 


2:3a2/oi)(l - 


2:702/01)' 




1 




I/0102 




(1 


— 2:4/0102) 


(1- 


2:4/ai02)(l - 


2:3/0102) 
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Thus applying the operator 61^^62.663 jSi,s to both sides of ([33]) gives 

'5l,5<52,6^3,7'54,8-F52 {xi,X2, X3, X4) 
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(1 — a;iaia2)(l — X2ai/a2)(l — a;3a2/ai)(l — 2:4/0102) 



(1 - X5aia2)(l - X6ai/a2)il - xra2/ai){l - Xslaia2) 



(34) 



Now we can easily recognize that ([M)) is precisely the constant term that MacMahon partition analysis 
would yield for the system SS2- This proves ([5^ . □ 



Note that to obtain the equality in (|34p we have used the simple fact that the divided difference 
operator and the constant term operator do commute. This is the fundamental property which is at 
the root of the present algorithm. This example should make it evident to have the following more 
general result (with double modified). 

Theorem 4.1. // Fs{xi, X2-, ■ ■ ■ , Xn) is the complete generating function of the Diophantine system 

bii bri ■■■ bin 



SS = 



brl br2 



brl br2 

n of the ( 
bii bi2 

brl br2 





Pi 




Cl 




_Pn_ 







s = 



then the complete generating function of the doubling of S defined by 
bii bi2 ■ ■■ bin 



bin 

b, 





Pi 






_P2n_ 





Cl - - 612 &1,; 

Cr — brl — br2 — ■ ■ ■ — brr, 



is given by the rational function 

Fssixi,X2, ■ ■ ■ ,X2n) = ^l,ri + l<52,n+2 ' ' ' <5n,2n^5 (a^l , 3^2 , ■ • ■ ,Xn)- 

This result combined with the next simple observation yields a powerful algorithm for computing 
a variety of complete generating functions. 

Theorem 4.2. Let Fs[xi,X2, ■ ■ ■ ,a;„) be the complete generating function of a Diophantine system 
S then the complete generating function Fss{xi,X2, . . . ,Xn) of the system SE obtained by adding the 
equation 

■S" = II ripi + r2P2 H h rnPn = s 

to S is obtained by taking the constant term 

Fse{xi,X2, . .. ,Xn) = a~''Fs{a''^xi,a'"'x2, ■ ■ ■ , a''"a;„) 



Proof. By assumption 



Now we have 



Fs{xi,X2, . . . ,a;„) = ^ : 



.P2 . . . „Pn 



pes 



a-'Fs{a'-^xi,a''^X2, . . . , a''"x„) = V • • • 2;^"a'-ifi+'-^P^+-+'^"f"-'* 

a*^ — ^ 



.P^ . . . -T.P" 



pess 

= Fssixi,X2, ■ . . ,a;„). 



□ 



These two results provide us with algorithms for (at least in principle) computing all the Hdd series 
Gk{q) as well as the Sdd series Wk{q)- 



14 A. GARSIA, G. MUSIKER, N. WALLACH, G. XIN 

Algorithm 4.3 (Hdd Case). ai) Initially compute the complete generating function for the Hdd 
problem for A; = 1. That is. compute the constant term 

Fi{xi,X2) = J- T^T ^ 

(1 - a;ia)(l - X2/a) 

ak) With Fk-i{xi, . . . ,X2k-i) from step bk-i), compute by divided difference 

FFk^l{xi, . . . ,X2k) = (5i,i+2''-i • • • '^2'=-i.2''^fc-l(a;i, • ■ • ,X2k-i)- 

bk) With FFk-i{xi, . . . ,X2k-i) from step ak), compute the complete generating function for the 
Sdd problem for k by the following constant term: 

Fk{xi,X2, ■ . ■ ,X2k) ^ FFk-i[axi,ax2, ■ ■ . ,ax2k-i,X2k-i+i/a, . . .,X2k/a) . 

This sequence of steps in Algorithm 14.31 can be terminated by replacing step bk) by 
bj^) The g-gencrating function Gk{q) is given by the constant term 

Gk{q) = FsSk-iiaq,aq,---,a.q,q/a,.-.,q/a) 

a" 

The steps up to bs) can be carried out by hand. For further steps we need a computer, and to 
carry out step bs) by computer we have to introduce one more tool as we shall see. Unfortunately 
Step be) appears beyond reach at the moment. 

It will be instructive to see what the first several steps give. 

ai) ^ 

(2^1, ,2:2) = 



a2) 



b2) 



FsSi{xi,X2,X3,X4:) = 

Fs^ (a;i,a;2,a;3,a;4) 



1 — X1X2 

(1 - XiX2X3Xi) 



(1 - 2;ia;2)(l - a;22;3)(l - a;ia;4)(l - xsx^) ' 
(1 - X1X2X3X4) 



(1 - a'^xiX2)il - X2X3){1 - a::ia:;4)(l - x^x^/a'^) 
1 



aa) 

Fss2 (x) 



(1 - X2a;3)(l - 2:1x4) 

(1 - 0:1x4x5X8) (1 - X2X3X6X7) 



(1 - XiX8)(l - X2X7)(1 - X3X6)(1 - X4X5)(1 - XiX4)(l - X2X3)(1 - X6X7)(1 - X5X8) 

ba) 

(1 - XiX4X5X8)(l - X2X3X6X7) 



Fs (x) — 

^ (1 - XiX8)(l - X2X7)(1 - X3X6)(1 - X4X5) 

^ (35) 



(1 — a^xiX4)(l — a^X2X3)(l — X6X7/a^)(l — xsXs/a^) 

We can compute this constant term in many ways. In particular we could use one of the MacMahon 
identities given by Andrews in [1]. But it is interesting to point out that our divided difference 
algorithm has already provided us (in step a2)) a formula we can use in step bg). In fact, the output 
of step 32) 

^ , s (I-X1X2X3X4) 

^55i(a:i,X2,X3,X4) = -Tj r-Tj ^ r-Tj r 

(1 - XiX2)(l - X2X3)(1 - XiX4)(l - X3X4) 

is the complete generating function of the system pi — p2 + Ps ~ P4 = 0, so by MacMahon partition 
analysis we should also have 

^5Si(2;i,a;2,X3,X4) 



(1 — axi)(l — X2/a)(l — ax3)(l — Xi/a) 
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This impUes that 



1 



(1 — a'^XiX4){l — a^a;2a;3)(l — X6XY/a'^){l — x^x^/a^) 



(1 - XiX2XzX4) 



(1 - XiX2){l - a:2a;3)(l - XiXi){l - X2.X4) 



(1 - XiX2X3XAXzX(iX7Xs) 



Xi- 


-fXiXi 


X3- 


^X2X3 


X2- 


->xeX7 


X4- 


■^X5Xs 



(1 - a;ia;4a;6X7)(l - X6a;7X2a;3)(l - a;ia;4X5X8)(l - X2X3X5X8) 

Using this in ([55)1 gives 

(1 - XiX4X5Xs){l - X2X3XeXY) 



Fssixi, ■■■,Xs) 



(1 - a;ia;8)(l - X2X7){1 - X3a;6)(l - ^4X5) 

(1 — a;iX2a;3X4a;5a:6a^7a;8) 



(1 — a;iX4X6a;7)(l — X6X7a;2a;3)(l — a::iX4X5X8)(l — a;2a;3a;5a;8) 
1 — X1X2X3X4X5XGXJXS 



(1 - xixs){l - X2a;7)(l - a::3a:;6)(l - a;4a;5)(l - a;ia;4a;6a;7)(l - a;2a;3a;5a;8) 
Replacing all the Xi by the single variable q, we thus obtain that 

Gi[q) = ^ , G2{q) = 2T2 ' <^3('7) 



l-g2' -^v^^ (l-q2)2' ^-^y-^J (1 _g2)4(i _g4)2 (1 - q^^ {1 - q^) ' 

Using the computer to carry out step gives 

1 + g2 + 21?-* + 36g6 + 7^ + 86gi° + 7^ + 36qU + 21q^^ + q^^ + q^" 

We shall see later what else has to be done to obtain G^{q). 

Our divided difference algorithm can also be adapted to compute the first 4 Sdd series as well. 
In fact, again due to the fact that divided difference operators commute with the constant term 
operators, we can also show that all the complete Sdd series can (in principle) be obtained by the 
following algorithm. 

Algorithm 4.4 (Sdd Case). ai) Initially compute the complete generating function for the Sdd 
problem for fc = 1. That is, compute the constant term 

1 - a2 

Wi{xi,X2) = -j r-j -r-r 

(1 - xia){l - X2/a) 

ak) With Wk-i{xi, . . . , X2k-i) from step bk_i), compute by divided difference 

WWk-l{xi, ...,X2k)= (5i,i+2'=-i • ■ ■S2k-l^2''Wk~lixi, ■ . .,X2k-l). 

bk) With WWk~i{xi, . . . , X2k-i) from step ak), compute the complete generating function for the 
Sdd problem for k by the following constant term: 

Wk{xi, X2, . . . , X2k) = WWk-i{axi, ax2, ■ ■ ■ , ax2k-i, X2k-i^i/a, . . . , X2k /a){l — a^) 

Note that similarly as for the Hdd-case, the sequence of steps in Algorithm 14.41 can be terminated 
by replacing step bk) by 

bj^) To obtain the generating function Wk{q) compute the constant term 
Wk{q) = WWk-i{aq, aq, . . . , aq, q/a, q/a){l - a^) 
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Only steps ai) and 3.2) can be carried out by hand. Though steps 3 and 4 are routme they are too 
messy to do by hand. But step 5 again needs further tricks to be carried out by computer. Step 6 
appears beyond reach at the moment. 

It wiU be instructive to see what some of these steps give. 



aij 

a2) 

WWiixi, 



Wi{xi,X2) 



1-4 

1 — XIX2 



,X4) 



1 — X2 — X2X4 — x1+ X1X2X4 + x'^XzXi — XiX2X3Xi + XiX2x1 + X2X^x'\ — Xix'^X^x'^ 
(1 - XiX2){l - X'iX2){l - XiXi){l - 0:3X4) 



W2{xi,X2,x:i,Xi) = 



This gives 

aa) 

WW2{xi,...,Xs) = 

W^ixi, ...,xs) 



W2iq) = 



1 - X2Xi - X3X4, + xj 
(1 - XlX4)(l - X2X3) 

1 



(large numerator) 



(1 - .Ti2;4)(l - XiXs){l - X2X3){1 - X2X7){1 - X3X6){1 - X4a;5)(l - X5Xs)il - XeXr) 

[large numerator) 



(1 - a;ia;8)(l - X2X7){1 - a;3a;6)(l - a;4a;5)(l - xiXiXeX7){l - a;2a;3a;5X8) 

ba) Notwithstanding the complexity of the previous results it turns out that to obtain 1^3(5) we 
need only compute the constant term 

1 l-a2 



Wsiq) 



(36) 



{l-q^a^){l-qya'^) 

To this end we start by determining the coefficients A and B in the partial fraction dccom 
position 

(1 -a2)a2 



1 



A 



B 



{\ — q^ 0^)10^ — (p) q^ \ — (f-a^ a^ — (p 



obtaining 



A^ 



B 



(l-a2)a2 



(a2 _ ^2) 
(1 -a2)a2 



(1 - l/<?2)/q2 



(l_q4) (1 + ^2)' 



(X-q^a"^) 

(the exact value of B is not needed) and we can write 
l-a^ 1 1 1 



q^ /a^ 



{l-q^a^){l-q^la^) q^ q'il + q') {\-a^q^) {I + q^) \ - q^ j a^' 
Thus taking constant terms gives 



1 



1 



(l-g2a2)(l_^2/„2) 

Using this in ([36|l we finally obtain 



,0 q^ 9^(1+92) 



= 



1 



1 



£14 



WW^(x\^X2t ■ ■ ^x\q) ~ (too large for typesetting) 
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h'^) Notwithstanding the complexity of the previous result it turns out that to obtain W4{q) we 
need only compute the constant term 



1 



(1 - q2)(l _ g4)2 (1 _ „2q4)(l _ ^4/^2)(l _ ^4^4)^^ _ <^4/„4) " 

To illustrate the power and flexibility of the partial fraction algorithm we will carry this out by 
hand. The reader is referred to [3] for a brief tutorial on the use of this algorithm. In the next few 
lines we will strictly adhere to the notation and terminology given in [3]. 

To begin we note that we need only calculate the constant term 



since we can write 



Now we have 



(1 - ax){l - x/a){l - a'^x){l - x/a'^) 
(l + g4)(l + 56) 



(l-g2)(l_g4)2 



X dq'). 



(37) 
(38) 



1 



x/a'^ 



{l-a^x){l-x/a^) {1 - a^x)ia^ - x) 
Thus (j37[) may be rewritten in the form 



1 



I — x'^ 1 — x/a^ 



C{x) 



1 



1 



(1-a) 



1 



(1-a) 



/ 2 

x/a 



(1 — ax){l — x/a) 1 — x/a^ 



(39) 



(1 — ax){l — x/a) 1 — a^x 

Note that in the first constant term we have only one dually contributing term and on the second we 
have only one contributing term. This gives 



(1-a) 



1 



(1 — aa;)(l — x/a) 1 — a^x 
(1 — a) x/c? 



(1-a) 



1 



(1 — aa;)(l — x/a) 1 — x/o? 
Using (gOl) and (gl]) in ([39]) we get 



(1 — ax) 1 — a^x 
(1 — a) x/c? 



{\-x) 



1 



(1 — x/a) 1 — x/ a^ 



(1 - .t2) 1 - x3 
-(1-x) a;2 



a— 1 / X 



(1 - a;2) 1 



(40) 
(41) 



C{x) 



(1-x) 



(1 



l-a;2 V(l-a;2) 1 



Together with ([35]), we get 



(1 -X2) 1-X^ 

l-g4 



(1-X2)(1-X3)' 



(l-g2)(l_q4)2 (l_q8)(l_^12) (I_52)(l_q4)2(l_g6)- 

We will see in section 4 what needs to be done to carry out step bg on the computer. 

The identities for W2{q),W3{q),W4{q) in ^ have also been derived in [5] by symmetric function 
methods from the relation In fact, all three results in ^ are immediate consequences of the 
following deeper symmetric function identity, (for a proof see [31 Section 2].) 



Theorem 4.5. 



Sd 



\h2d 



where EO4 denotes the set 0/ partitions 0/ length 4 whose parts are > and all even or all odd. 

Note that the Kronecker product identity 

{sd.dQ Sd,dQ SdAQ Sd,dQ Sd,d , S2d) = {sd.dQ Sd^dQ Sd,d , Sd,dQ Sd,d)- 

suggests obtaining W5{q) by means of a combinatorial interpretation of the coefficients of the Schur 
function expansion of the Kronecker product Sd,d © Sd,d © Sd,d- However, to this date no formula has 
been given for these coefficients, combinatorial or otherwise. 
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5. Solving the Hdd problem for fc = 5 

This section is divided into four parts. In tlic first subsection we start with our computer findings 
and end by giving a combinatorial decomposition that works nicely to obtain i^3(a;). In the second 
subsection, this decomposition is described algebraically and, together with group actions, turned into 
manipulatory gyrations that will be used to extract G^{q) and W^{q) out of our computers. In the 
third subsection, by combining the idea of decomposition and the method of divided difference in 
Section m we give our best way that reduce the computation time for G^[q) and W^{q) down to a few 
minutes. In the final subsection, we give our first algorithm to obtain G^{q) and ^5(9). 

5.1. A combinatorial decomposition for ^3(2;). Our initial efforts at solving the Hdd an Sdd 

problems were entirely carried out by computer experimentation. After obtaining quite easily the 
series G2(g), 6*3(9), G/^{q) and W2{q), Wz{q)^ Wi{q), all the computer packages available to us failed 
to directly deliver G^[q) and W^{q). 

The computer data obtained for the Hdd problem for fc = 2, 3 were combinatorially so revealing 
that we have been left with a strong impression that this problem should have a very beautiful 
combinatorial general solution. Only time will tell if this will ever be the case. To stimulate further 
research we will begin by reviewing our initial computer and manual combinatorial findings. 

Recall that we denoted by Td the collection of all d-subsets of the 2d element set Vl2d- We also 
showed (in Theorem 13. ip that the coefficient md{k) in the scries Gk{q) = J2d<o1^'^^d.ik) counts 
the number of orbits under the action of the symmetric group S2d on the fc-fold cartesian product 
X J-'d X ■ ■ ■ X Td- Denoting by {Ai, A2, . . . , A^) a generic element of this cartesian product, then 
each orbit is uniquely determined by the 2'' cardinalities 

where for each 1 < i < fc we set 

^■'={ 4 if^;^?; ^^"'^ '=A,=n2d/A,). 

It is also convenient to set A^-^^^^^... = A^^ n n • • • n A'j!' . This given we have seen that the 
condition {Ai, A2, . . . , Ak) € is equivalent to the Diophantine system 

2ei)Pei,e2,--,ffc = 0, 
2efc)Pei,e2,-,efc = 0, 

together with the condition |f22d| = 2d, that is J2ei=o J2e2=o ' ' ' J2lk=oP'^i,'^2,--- ,>^k = 2d. 

There are several algorithms available to solve such a system. See for instance [3 Chapter 4.6]. 
The algorithm we used for our computer experimentations is the MacMahon algorithm which has 
been recently implemented in MATHEMATICA by Andrews, Paule and Riese and in MAPLE by Xin 
using the partial fraction method of computing constant terms. 

The former can be downloaded from the web site 

h ttp:/ / www. rise, uni-linz.ac. a t/ research / combina t/ software/ Omega / 

and the latter from the web site 

^ h ttp:/ / www, combina torics.net. cn / homepage/xin /maple/ el 12. ra^ 

For computer implementation we found it more convenient to use the alternate notation adopted in 
Remark [131 That is 

Sk = II PiVi +P2V2 + ■■■ + P2^V2k = . (42) 

These algorithms may yield quite a bit more than the number of solutions of such a system. For 
instance, in our case letting Ck denote the collection of solutions of the system Sk, the "Omega 
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package" of Andrews, Paule and Riese should, in principle, yield the formal power series 

Fk{xi,X2,...,X2k) = ^ X^^ X^^ ■ ■ ■ . 

(pi,P2,---,P2fc)eCfc 

It follows from the general theory of Diophantine systems that Fk{xi, X2t ■ ■ , X2k) is always the Taylor 
series of a rational function. 

Now for ^2 and S^, the Omega package gives 

F2{xi,X2,x-i,x^) = (43) 

(1 - XiXi)[l - X2X'i) 

I - X2X3Xr,Xs,XlXiXQX^ 

Fz{xi,X2,...,x^)^ — — — — — -. (44) 

(1 - a;ia;8)(l - X2X7)(l - a;3X6)(l - X4^x^)(l - X2XzX5Xs,){l - xiXiX^xj) 

But this is as far as this package went in our computers. However we could go further by giving up 
full information about the solutions and only ask for the series 



Gk{.q) = Fk{xi,X2, ■ ■ ■,X2k] 



which can be computed from its constant term representation in ([8]). For example, the program Latte 
by De Loera, Hcmmeckc, Tauzcr, Yoshida, which is available at 

http:/ / www. math. ucdavis. edu /~ latte/ 

computed the Gi{q) series in approximately 30 seconds. However, this is as far as Latte went on our 
machines. We should also mention that all the series Gk{q) and Wk{q) for A: < 4 can be obtained in 
only a few seconds, from the software of Xin by computing the corresponding constant terms in ([5|) 
and dl]). 

To get our computers to deliver Gz{q) and W^{q) in a matter of minutes a divide and conquer 
strategy had to be adopted. More precisely, these rational functions were obtained by decomposing 
the constant terms ([8]) and ([9]) as sums of constant terms. This decomposition had its origin from 
an effort to find a human proof of the identities in (j43|) and (|44p . More importantly, the surprising 
simplicity of (|43p and (|44p required a combinatorial explanation. Our findings there provided the 
combinatorial tools that were used in our early computations of G^{q) and Wz{q). This given, before 
describing our work on these series, we will show how to obtain (j43p and (|44p entirely by hand. 

Let us start by sketching the idea for k = 2. Beginning with 

^ ^ Pi + P2 - P3 - P4 = 
^ Pi - P2 + P3 - Pa = 

we immediately notice that (1, 0, 0, 1) and (0, 1, 1, 0) are solutions. Set 

a = min(pi,p4) and b ~ min(p2,P3)- 

It is clear that the following difference must also be a solution. 

(91,92,93,94) = (P1,P2,P3,P4) - {a,b,b,a) = (pi - a,p2 - b,p^ - b,p4 - a). 

Now (71(74 = and (7293 = 0. This gives us four possibilities for ((71, 52, 93, 94): 

(0,0,x,2/), iO,x,0,y), (x, 0,27,0), (a:, 27,0,0), (45) 

for some nonnegative integers x,y. Testing the first equation of ^2 immediately forces the first and 
last in ([^5]) to identically vanish. Similarly, the second equation of 52 yields that the second and 
third in P5|) must also identically vanish. This proves that the general solution of 1S2 is of the form 
(a, b, b, a). We thus rcobtain the full generating function of solutions of 52) 

i^2(a;i,a;2,a;3,a;4) = y^^^i^^xl^l = -!- -. 

(1 — a;iX4)(l — X2a;3) 

a>0 6>0 ^ i *M ^ a/ 

It turns out that we can deal with S3 in a similar manner. Again we begin by noticing the four 
symmetric solutions 

(1,0,0,0,0,0,0,1), (0,1,0,0,0,0,1,0), (0,0,1,0,0,1,0,0), (0,0,0,1,1,0,0,0). 
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Next we set 

a = min(pi,p8), b = inm{p2,P7), c ^ mm{p3 , pe) , d = min(p4,p5), 
and by subtraction we get a solution 

(91,92,93,94,95,96,97,98) = (Pi,P2,P3,P4,P5,P6,P7,P8) ^ {a, b, c, d, d, c, b, a) (46) 
with the property qiqg-i = for 1 < i < 4. It will be good here and after to call the set 

{i e : Pi>l} 

the support of the composition (pi,p2, ■ ■ ■ ,Pn)- This given, we derive that the resulting composition 
in ([46]) will necessarily have its support contained in at least one of the following 16 patterns. 

(0,0,0,0,*,*,*,*), (0,0,0,*, 0,*,*,*), (0,0,*,0,*,0,*,*), (0,0,*,*, 0,0,*,*), 
(0,*,0,0,*,*,0,*), (0,*,0,*,0,*,0,*), (0,*,*,0,*,0,0,*), (0,*,*,*, 0,0,0,*), 
(*, 0,0,0,*, *,*,0), (*,0,0,*,0,*,*,0), (*,0,*,0,*,0,*,0), (*,0,*,*,0,0,*,0), ^ '> 
(*,*,0,0,*,*,0,0) , (*,*,0,*,0,*,0,0) , (*,*,*,0,*,0,0,0) , (*,*,*,*, 0,0, 0,0). 

Unlike the case k = 2 not all of these patterns force a trivial solution. To find out which it is helpful 
to resort to a Venn diagram imagery. To this end recall that a solution of 1S3 gives the cardinalities of 
the 8 regions of the Venn diagram of three d-subsets Ai,A2, of (see Figure [5]). 




Figure 2 . The Ven diagram for . 

In Figure [31 each pattern is represented by a Venn diagram where in each region A\^ n A^^ n A^^ that 
corresponds to a * in the pattern we placed a black dot. That means that only the regions with a 
dot may have > cardinality. The miracle is that all but the two patterns (0, *, *, 0, *, 0, 0, *) and 
(*, 0, 0, *, 0, *, *, 0) can be quickly excluded by a reasoning that only uses the positions of the dots in 
the Venn diagram. In fact, in each of the excluded cases, we show that it is impossible to replace 
the dots by > integers in such a manner that the three sets ^1,^2,^43 and their complements 
'^Ai, '''A2, ''A3 end up having the same cardinality (except for all empty sets). 




Figure 3. The 16 support patterns for ^3. 



The reasoning is so cute that we are compelled to present it here in full. In what follows the j 
diagram in the i"^ row will be referred to as "Dij": 

(1) Dii, Di4, DiQ, D23, D25, and Z?28 can be immediately excluded because one of Ai, A2, A3, 
Al, A2 or A3 would be empty. 
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(2) In Di5 the dot next to 8 should give the cardinahty of (say d) and then the dot next to 
the 2 should also give d. But that forces the dots next to 5 and 6 to be 0, leaving ^3 empty, 
a contradiction. The same reasoning applies to -D12, £'13, Dis: D21, D^i, £'267 ^-nd ^27- 

That leaves only the two diagrams Dn and D22 which clearly correspond to the two above mentioned 
patterns. Now we see that for D22 we must have the equalities pi + p4 = pi + pe = pi + pr = 
Pe +P7- This forces pi = Pi = pe = Pi- In summary this pattern can only support the composition 
(m, 0, 0, u, 0, u, u, 0). The same reasoning yields that the diagram Dn can only support the composition 
(0, V, V, 0, V, 0, 0, v). It follows that the general solution of ^3 must be of the form (a, 6, c, d, d, c, b, a) + 
{u, V, V, U, V, u, u, v). 

Now recall that after the subtraction of a symmetric solution we arc left with an asymmetric 
solution. Thus to avoid over counting we must impose the condition uv = Q. This leaves only three 
possibilities m = w = 0, u > 0, w = or w = 0, w > 0. Thus 

{xiXiXQXiy- + ^'^{x2X^XzXf,)'"^ 
a>0 6>0 c>0 d>0 u>l v>l 

_ 1 A ^ XiXjXQXl ^ X2X3X5XS 

(1 - a;ia;8)(l - X2X7)(1 - X3a;6)(l - a;4a;5) \ 1 - a:iX4X6a;7 1 - X2X3X5XS 
which is only another way of writing (|44p . 



5.2. Algebraic decompositions and group actions. It is easy to sec that the decomposition of a 
solution into a sum of a symmetric plus an asymmetric solution can be carried out for general k. In 
fact, note that if < i < 2'° — 1 has binary digits £462 •••£/£ then the binary digits of 2*^ — 1 — i are 
£162 ■ • ■ Efc (with e = 1 — e). Thus we see from (|42|) that in each equation pi and p2'=+i-i appear with 
opposite signs. This shows that for each k > 2 the system Sk has 2*^"^ symmetric solutions, which 
may be symbolically represented by the monomials XiXii for i = 1, . . . , 2*^"^, where we use (and will 
often use) i' to denote 2*^ + 1 — z when k is fixed. 

Proceeding as we did for 52 and 53 we arrive at a unique decomposition of each solution of Sk into 

(pi,P2, • ■ • ,P2>') = {U1,U2, . . ■,U2,Ui) + iqi,q2, ■ ■ • ,920 

with the first summand symmetric and the second asymmetric, that is Ui = Ui' and QiQi' = for 
1 < i < 2''^^, and thereby obtain a factorization of Fk{x) in the form 

Fkix)=(ll—^ (48) 

\ ^ XiXi' J 

with F^{x) denoting the complete generating function of the asymmetric solutions. 



This given it is tempting to try to apply, in the general case, the same process we used for fc = 3 
and obtain the rational function F^{x) by selecting the patterns that do contain the support of an 
asymmetric solution. Note that the total number of asymmetric patterns to be examined is 2^ 
which is already 256 for fc = 4. For fc = 5 the number grows to 65, 536 and doing this by hand is 
out of the question. Moreover, it is easy to see, by going through a few cases, that even for fc = 4 
the geometry of the Venn Diagrams is so intricate that the only way that we can find out if a given 
pattern contains the support of a solution is to solve the corresponding reduced system. 

Nevertheless, using some inherent symmetries of the problem, the complexity of the task can be 
substantially reduced to permit the construction of ^5(5) by computer. To describe how this was 
done we need some notation. We will start with the complete generating function of the system 
as given in Remark l2.5i that is 



Fk{xi,X2, ■■■,X2k) =Y\_ 
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riiLi '^i with £162 •••efe being the binary digits of z — 1. Note that since (as we 



where Ai 

previously observed) the binary digits of 2^^ 
that 

1 •^i^i' 



I are £162 • 
1 

1 ^ XiAi 



Thus combining the factors containing Ai and Aii we may rewrite 

Fkixi,x2,...,x2k) = n T^z— n 



■ Efe , we have A 

Xi'/A 



l/Ai. It then follows 



Xi' I Ai 

in the form 



1 



Xi' I Ai 



{\-XiAi) (l-xvjA,) 



(49) 



aVa^'-'-aV. 



i—\ i—\ 

Comparing with (|48p we derive that the complete generating function of the asymmetric solutions is 
given by the following sum. 



SC[1.2''-i 



where 



Fs{x) 



(1-.T,A,) 



I Ai 



Xi' I Ai) 



(50) 



(51) 



In this way we have described our decomposition algebraically. Using notation as of ([55]) , we can see 
that Fs{x) is none other than the complete generating function of the reduced system 



PiVi + ^Pi'Vi' 







with the added condition that pi/ > 1 for all i E S. 

Note that for fc = 3 the summands in ([50]) correspond precisely to the 16 patterns in (|47|) with the 
added condition that the in position i > 5 should represent > 1 in the corresponding solution 
vector. This extra condition is precisely what is needed to eliminate overcounting. 

Perhaps all this is best understood with an example. For instance for fc = 3 the patterns 
(*,0,0,*,0,*,*,0) and (0, *, 0, *, 0, 0, *) 
were the only ones that supported an asymmetric solution represent the two reduced systems 



{23} 



^{14}= Pi - P4+P6 - P7 = S 
Pi - Pi- P6 +P7 = 

and correspond to the following two summands of (jSOp for fc = 3 

1 1 Xea2/aia3 0:703/0102 



P2+P3~P5- 
P2-P3+P5- 
-P2 +P3+P5 



PS = 
P8 = 
~P8 



^{1,4} (2;) 



Fi 



1 — Xiaia2a3 1 — a:;4ai/a2a3 1 — X6fl2/aia3 1 — 2:703/0102 
1 1 2:50203/01 2:8/010203 



'{2,3} (2:) — -, — 1 — -; 1 — -; -. 

1-2:20102/031-2:30103/021-2:50203/011-2:8/010203 



2:1X42:62:7 
1 — 2:12:42:52:7 

2:2X32:52:8 
1 - 2:22:3X52:8 



(52) 
(53) 



A close look at these two expressions should reveal the key ingredient that needs to be added to 
our algorithms that will permit reaching fc = 5 in the Hdd and Sdd problems. Indeed we see that 
i^li 4}(a:) goes onto -F'{2.3}(2:) if we act on the vector (xi,X2, • • ■ ,2:8) by the permutation 

'1 2 3 4 5 6 7 
,3 4 1 2 7 8 5 6, 



(54) 



and on the triple (01,02,03) by the operation 02— s-Oj^^. In fact, a is none other than an image of the 
map (ei, £2, e3)^(ei, £2, £3) on the binary digits of 0, 1, . . . , 7, as we can easily see when we replace 
each i in ([5^ by the binary digits of i — 1 

000 001 010 Oil 100 101 110 111 
010 Oil 000 001 110 111 100 101 
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What goes on is quite simple. Recall that solutions p of our system Sk can also be viewed as 
assignments of weights to the vertices of the /c-hypercubc giving all hyperfaces equal weight. Then 
clearly any rotation or reflection of the hypcrcube will carry this assignment onto an assignment with 
the same property. Thus the Hypcroctahedral group will act on all the constructs we used to solve 
Sk- 

To make precise the action of Bk on [1, 2'^] we need some conventions. 

(1) We will view the elements of Bk as pairs (a, 77) with a permutation a = (ai, Q!2, . . . afe) G Sk 
and a binary vector ry — (rji, ?/2, ■ • ■ , 'yfc)- 

(2) Next, for any binary vector e = (ei, 62, . . . , e^) let us set 

(a, r])e = (cai + + '72, • ' ' , + Vk) (55) 

with "mod 2 " addition. 

(3) This given, to each element g = (a, if) G Bk there corresponds a permutation a{g) by setting 

1 2 ••• 2^ 

,(Ti (72 ••• 0-21- 



<9) 



where cr, = j if and only if the A;- vector e = (ei, 62, . . . , e^) giving the binary digits of z — 1 is 
sent by (|55p onto the /c- vector giving the binary digits of j — 1 . In particular we will set 

g(xi,x2, ...,X2k) = {xa,,Xa2,- ■ -^x^^J. (56) 
(4) In the same vein we will make Bk act on the fc-tuple (ai, a2, . . . , Cfc) by setting, again for 
9 = {a, 77) 

g{a,,a2,...,ak) = (aj,7'"\ ai;'"% . . . , a^;''"=). (57) 
With these conventions we can easily derive that gXiAi = Xa- A„. . Thus 

TT 1 _ TT 1 TT ]_ 

^^l-XiAj nn n ^^l~X„,A„, nn n ^ ^ 1 - Xi 



iAj 



a';a);---a", 



from which we again derive the Bk invariance of the complete generating function Fk{xi, X2, ■ ■ ■ , X2k). 

If we let Bk-i not only act on the indices 1, 2, ... , 2^^"^, but also on 1', 2', . ■ • , 2*^"^' by ai' = <t[. 
Then Bk-i permutes the summands in ([5D|) as well as the factors in the product 

1 XiXi' 

'l—l 

Note further that if we only want the g'-series Gk {q) we can reduce ([501) to 

Gt{q)^ E Gs{ql (58) 
SC[l,2''-i] 

where Gs{q) ~ Fs{x)\^^^. But if for some g € Bk-i we have 

(Xo-i , , . ■ . , Xcr^ J = Fs:,ixi,X2, ■ ■ ■ ,X2k), 

then replacing each Xi by q converts this to the equality Gsi (q) = (q) ■ That means that we need 
only compute the constant terms in (|58[) for orbit representatives, then replace (j58[) by a sum over 
orbit representatives multiplied by orbit sizes. More precisely we get 

G^(g) = 5]m.Gs.(g), (59) 

i 

where the sum ranges over all orbits and rrii denotes the cardinality of the orbit of the representative 
Fsi{x). In the computer implementation we obtain orbit representatives as well as orbit sizes, by 
acting with Bk-i on the monomials Ms = Ilies^Ji. 

Thus for = 3 we found that the 16 summands in (|50p break up into 6 orbits but only 2 of them do 
contribute to F^. They corresponds to the monomials 1 and X1X4 with respective orbit sizes 1 and 2. 
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The orbit representative that corresponds to 1 is simply the case S = 4> in ((5T|l and that corresponds 
to X1X4 is given in ([5^ . 

Thus from dSH), and (gl]) we derive (again) that 
^3(9) = n 3T4 1 



(l-q2)4V l-qi) (l-g2)4i_g4- 

For fc = 4 we have 2^ = 256 summands in ([50|l with 22 orbits but only 11 of these orbits do contribute 
to F^. The number of denominator factors for each term is 8 which is still a reasonable number for 
the partial fraction algorithm. The formula for ^4(0;) obtained this way can be typed within a page, 
but we would like to introduce a nicer F^ix) using the full group Bk instead of Bk^i, as we will do 
in the next paragraph. For fc = 5 we have 2^^ summands in ([50)1 with 402 orbits but only 341 orbits 
do contribute to F^. The number of denominator factors for each term is 16 which is out of reach 
for the partial fraction algorithm to obtain F^{x). Nevertheless, in this manner we can still produce 
G5{q) in about 15 minutes. 



The decomposition in ((58|) is only Bk-i invariant, and it is natural from the geometry of the 
hypercube labclings, to ask of a Bk invariant decomposition. To obtain such a decomposition of F^^x) 
we will pair off the factors containing Ai and Aii by means of the more symmetric identity 

1 XiX2>'-\-i—i f ^ , x^Aj^ ^ XiiAii 



{l~ x,A,){l - x2k+i_.,A2k+i_i) \ l~x^A, l-XifAe 
and derive that 

Fk{x) = Ps,t{x), 

SUTC[1, 2*^-1] 



XiAi \ / T— r x'j^/Ai 



where S and T are disjoint and 

Note that every pair (5',T) should be identified with the set S U {i' : i e T} C [1,2'=] when applying 
the action of Bk- 

Example 5.1. For fc = 3 we have 3^ = 81 summands with 9 orbits but only 2 orbits do contribute 
to F^. The two orbits corresponds to the monomials 1 and xiX4XqX7 with respective orbit sizes 1 
and 2. The orbit representative that corresponds to 1 is simply the case F^^^ = llaOo^ai^a^ ^ 1 ^^^d 
that corresponds to XlX4Xex^ is 



xiAi X4A4 xqAq xjA 



7 



F{l,4},{2,3}ix) ^ _ ^^^^ ^ _ ^^^^ ^ _ ^^^^ ^ _ ^^^^ 

Therefore, we can rcobtain G3{q) by using ([59]) as follows. 



XiX4XqXj 
1 — XiX4XqXj 



1 , 2g^ ^ l + g4 



GM = T. ?T4(1 



(l-g2)4^ l_g4/ (1 -q2)4(i _q4)- 

Example 5.2. For fc = 4 we have 3* = 6561 summands with 62 orbits but only 10 orbits do contribute 

(1) 1 

X\X\^X4X\4 



to F^ . We obtain the following complete generating functions for the 10 orbit representatives: 



(24) 



1 — X\X\^X4X\4 



XiqXt {xgfxeXj 
1 — XiQXTXg'^XQX4 

rnci\ a;i5a;3a;7 (2^12)^ {xgf {xef 

(1 - a::i2X7a;9a;6) (1 - xisXsXuXgXa"^) 
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(96) 



(192) 



(1 - a;2X7Xiia;i4) (1 - a;i6a;5X2a;ii) 

XgXioXi (Xi)'^ (xi^f (xs)^ {xuf 



(1 - a;ia;i5X4a;i4) (1 - X15X5X1QX4) (1 - a;i5X5a;9a;42a::i4) 



(64) 

(32) 



.g^s XeXl6X4 {xsf {x^f {Xi5f {Xiof 

(1 - a;i5X5XioX4) (1 - xiex^xaXio) (1 - a;i5a;3a;ioa;6) 

X3XrX4 {xefxi (xg)^ {Xl2f {xisf 

(1 - a;ia;i5a;i2a;6) (1 - xi2XrXgX6) (1 - xisXgXeXi) (1 - xi5X3Xi2a;9a;62) 

{xi^f ixi2)^ a;ia;3a;2a;6X7a;8 (l - a;ia:;2a;3a;8a;i2^a;7a;i3^a;6) 



(1 - xiX8a;i2Xi3) (1 - X2Xi2a;7a::i3) (1 - cc3a;i2a:;i3a::6) (1 - a;ia;i22x7a;i3a;6) (1 - a;2a;3a;8a:i2a;i3^) 

(8) a;4a;5a::3a;6a;9a;ioa;i5Xi6 (l - 2 xi^xwXzXsXioXgXeXi + xi^'^xie'^x^'^xs'^xio'^XQ'^XQ'^Xi'^) 
(1 - xiexsxgxe) (1 - xiex^xgXi) (1 - a;i5X9X6a;4) (1 - a;i5a;5Xioa;4) (1 - a;i6a;5X3a;io) (1 - xisXsXioXe) 

Here the numbers in parentheses give the respective orbit sizes. 
Replacing all the Xi by q and summing as in (|59p . we obtain 

l + q^ + 2lq'^ + 36g6 + 7V + 86gi° + 7Aq^^ + 3691-^ + 2lq^^ + + g^o 



Gi{q) 



(l_g2)7(l_^4)4(i_^6) 



We should mention that the partial fraction algorithm delivers this rational function in less than a sec- 
ond by directly computing the constant term in ([8]) for fc = 4. We computed the above representatives 
because it contains more information and can be used for an alternate path to Gs {q) . 

Computing the orbit representatives for k ~ 5 requires the construction of the 2^ x 5! — 3840 
elements of B5 and examining their action on the 3^^ = 43046721 symmetric supports. This took a 
few hours on our computers. We found in this manner that the 43046721 summands in (|58p break up 
into 15418 orbits and of these 6341 contribute to the sum. Most of the orbits have denominators of 
less than 16 factors. It also took about 15 minutes to persuade MAPLE to deliver 65(5) in the form 
displayed in the introduction. 

It turns out that the same orbit reduction idea can also be used to compute W5{q), but much more 
complicated. Let us explain the details in the next subsection. 

Remark 5.3. It is interesting to point out that computing complete generating functions for orbit 
representatives of summands in (j50p yielded as a byproduct orbit representatives of the extreme rays 
of our Diophantine cone for fc = 4 and fc = 5. Note that for fc = 3 the representatives can be directly 
derived from our hand computation, there are only two and the corresponding Venn Diagrams are 




Here the regions without numbers are empty. The number 1 indicates that the region has only one 
element. For fc = 4 we found that there are only three orbits, containing 24, 8 and 16 elements 
respectively, the corresponding diagrams are depicted below. 
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Note, for k ~ A each Venn diagram is depicted as a pair of Venn diagrams of fc = 3. The first 
member of the pair renders the Venn diagram of Ai H A2, Ai H A3, Ai A4 and the second member 
renders the Venn diagram of ^Ai n A2, ''Ai n ^3, "Ai n A4. 

For fc = 5 we found that there are 2712 extreme rays which break up into 9 orbits. We give in 
Figure [1] a set of representatives depicted as assignments of weights to the vertices of the 5 dimen- 
sional hypercube. We imagine that the vertices of this hypercube are indexed by the binary digits 
of 0, 1, 2, ... , 31 with 00000 the vertex at the origin and 11111 giving the coordinates of the opposite 
vertex. In Figurc[3]each hypercube is represented by two rows of two cubes. The cubes in the first row, 
from left to right, have the vertices labeled with the binary digits of 1 to 16 (minus 1) and the cubes 
in the second row have the vertices labeled with the binary digits of 17 to 32 (minus 1). The vertices 
here have possible weights 0, 1, 2, 3 and, correspondingly, are surrounded by 0, 1, 2, 3 concentric circles. 
The integer on the top of each diagram gives the size of the corresponding orbit. 



7\ 



V 



Figure 4. Representatives of extreme rays for /c = 5. 



Each of the corresponding solutions of our system is minimal, that is, it cannot be decomposed 
into a non-trivial sum of solutions. But we found that there are also 480 minimal solutions that do 
not come from extreme rays. The latter break up into two orbits, with representatives depicted in 
Figure [S] 

5.3. Our fastest way for G^lq) and Wz{q). With the notations in the previous subsection and 
Section [4] handy, we can describe our best way to obtain G5 (g) and W5 {q) . 
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V 



Figure 5. Representatives for minimal solutions but not extreme rays. 



Let us explain the idea for k = 5. In Example 15.21 we have obtained for F^{x) 10 orbit representa- 
tives with corresponding orbit sizes. Denote them by Ri{x) the representatives and nii the orbit sizes 
for i = 1, . . . , 10. From this we can give explicit formula of F^{x) and hence of ^4(0;) with the help 
of B4 action as follows. 



Ff{x) 



Ri 



geBi 



n«=i(l -2;»2;i7-i) 



Applying Algorithm 14.31 a.»;). b.»;) to ([50]) . we can obtain ^5(0;) by multilinearity. 



10 



Ri 



j = l,2,...,16 



£Cj=£Cja,£Cii3+j=xi6+j /a 



(60) 



(61) 



where we have used the straightforwardly checked fact: for any rational function R{xi, . . . ,xig) and 
g E Bk, it holds that 

^1,17 ■ ■ -^16, 32 gR{xi, ■ ■ ■ ,xie) = g5i,n ■ ■ ■ 5i6,32-R(a:^i, • ■ • ,a;i6), 
where g is extended to permute also indices 16 + j by ^(16 + j) = 16 + g{i) for j = 1, . . . , 16. 
Substituting Xj = q for all j into (j6ip gives 

10 / n \ J = l,2,...,16 



Ri 



Jl6,32:;— g 



4=1 



(62) 



That is to say, we only need representatives of Fk-i{x) together with orbit sizes to compute Fk{x), 
and this clearly extends for general k. Using (p^ . we can persuade Maple to deliver G5{q) as in (fTU)) 
in about 12 minutes. 

The orbit reduction idea for Gr:,{q) works in a similar way for W^^q). In fact, we can carry out almost 
verbatim the same steps that yielded the orbit decomposition of the complete generating function 
Fk{xi, X2, ■ ■ ■ , X2k ) to obtain the complete generating function Wfc(xi, a;2, . . . , X2k) as we shall define. 
Recall that the Wk{x) was originally defined in (|23|) as the constant term 



k 2" ^ 

Wk{xi,X2,...,X2k) = n(^^"i)n^i T 



j=i 1=1 

To carry out its decomposition we need only observe that if we let 



Wkixi,X2, ■ ■ ■,X2k) 



^n(i-s^)(i-ar ' 



^ n 1 _ ^^A, 



(63) 



(64) 



then 



Wk{q) = Wk{,q). 

The reason for this is that when all the Xi are replaced by we can easily show that the constant 
term in (|63p is not affected if we replace any ai by . Thus if we average out the right hand side of 
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((63)) over all these interchanges the result will be simply the right hand side of ((64l) due to the simple 
relation 

l^'^±f^ = \{l-aj){l-af). 

Now ([M)) brings to evidence that Wk{x) is Bk invariant while Wk{x) is not. Symmetrizing Wk{x) 
gives Wk{x). We can obtain either a Bk-i invariant decomposition or a Bk invariant decomposition 
of Wk{x) just as for Fk{x). 

The orbit reduction can also be used to considerably speed up steps at) and bj^) in Algorithm 14.41 
of the divided difference. The idea is similar as for the computation of G^{q), but is much harder to 
be carried out. 

To be clearer, we note that in step bj^) we do not need the complete generating function Wk-i{x). 
One way is to replace it by the more symmetric Wk-i{x). We have 

Wk-i{x) - _^ W?s,t(x), 

where S and T are disjoint as before. We only need to find orbit representatives 

with respective multiplicities mi, TO2, ■ • ■ , 'titv, since from them we can rebuilt Wk-i{x)^ just as in 
(|60)) . Then in step ak) we can replace Wk-i{x) by the sum 



1 ^ 
WU [x) = -—^ m.,Ws^,T. {x) 



and, with a similar reasoning as for Gk{q), obtain 

Wk{q) = y "^^Ul,l+2'=-l • • •'52'=-i.2'=W^Si,T.(2;) (1 - a^) 



(65) 



When working with (g) , we need an analogue of the collection of orbit representatives together 
with orbit sizes as in Example [521 Although Maple gives such a collection, we find it too complicated 
to be handled by Maple when using (p5)) . 

We find a way to avoid this problem. The idea is that in a formula like the Ri need not be 
chosen to have combinatorial meanings. This is best illustrated by the fc = 3 case. We can clearly 
see the advantage of orbit reduction in producing a compressed version of W-k{x). For fc = 3, the B^ 
decomposition will give 9 orbits with only 7 of them contributing to W^^x). We thus get 

~A 1 ■r-^ / TJraonomials 
"3 yx) = Tm / 9 9 monomials H 

03 7^ ' V 1 ~ XiXiX6X7 

The actual formula is a little complicated and its combinatorial meaning is not significant, but it 
is good enough for us to use the divided difference algorithm to compute W4{q). From this, by 
symmetrizing and re-choosing representatives, we obtain a simpler representative. Namely we end up 
obtaining that 

1 / 2 - 6x1X7 -xi"^ + 6x1X4^ xr-xi'^XA^xr^ 

S3 Tl^ \ (1 - X1XGX4X7) 



which can also be used in our divided difference algorithm. Originally we hoped that this formula 
would enable us to compute W^^q) entirely by hand, but we were unable to do so. 

For fc = 4, directly using the S4 decomposition gives us 62 orbits with 27 of them contributing to 
W4{x). The representatives obtained this way are too complex for further computation since several 
of them have thousands of monomials in their numerators. The similar idea of symmetrizing and 
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re-choosing applies to give us 10 reasonably simple representatives for W4{x), but typesetting them 
will take several pages. Nevertheless we are able to use them in the divided difference algorithm. 

Having noticed that for k = 2,3,4 the divided difference algorithm reduced the computation of 
Wk{q) to a rather simple constant term evaluation, we tried to see what it gave for fc = 5. Adding the 
contributions of these 10 representatives, before taking the constant term, yielded a rational function 
of the form 

1 

(1 - g2) (1 _ ^4)4 _ ^6) _ f^) (1 - a2g2) (i _ a^)' (1 _ a^qif 

357 monomials 

a^q^f (l - f^) (1 - (i - f^) (1 - aV) (l - (1 - a^<Z«) ' 

It turns out that this is actually a rational function in and a^. Replacing q by q^^^ and a by a^/^ 
and then taking constant term in a, we can obtain W5{q^^^). Using this approach Maple can deliver 
Wz{q) in only about 5 minutes in total which is the shortest time we have been able to compute this 
series. 



5.4. Our first algorithm to obtain G^lq) and W5{q). Before closing it will be worthwhile to include 
a description of the first algorithm that was used to obtain G5{q) and Wr^{q) since it contains another 
trick that clearly shows the flexibility afforded by the partial fraction algorithm in the computation 
of constant terms. 

In this approach we begin by replacing our system Sk by a system Sj^ which has the same cone of 
solutions. To describe the new system we will use the fc-tuple of sets model. The idea is that originally 
we got Sk by equating the cardinality of each set to the cardinality of its complement obtaining 





l^ll = 




Sk = 


\A2\ = 






\Ak\ = 


'Ak 



Now it is quite clear that this is equivalent to set 



S'k - 



\Ai\=d 

\A2\=d 

\Ak\ = d 
Mil 

For instance, using the binary digit indexing of the variables, for k 
system of 4 equations in 9 unknowns 



3 this results in the following 



Paoa 
Pooo 
Pooo 



Pool 
Pool 



Paw 



Pow 



Poll 



+ 



PlOO 
PlOO 
PlOO 



+ PlOl 



PlOl 



PllO 
PllO 



Pill 



-d 
-d 
-d 
- d 










This given, our rational function G3{q) 
constant term 

G3{q,t) = ^ ^ 



G3{q,l) may be also obtained by taking the following 



1 



1 



1 - 9010203 1 



qaia2 1 
1 



gaioa 1 
1 



qai 



1 



1 



1 



a'-;a'Aa'ia^: 



1 — 5020304 1 — (702O4 1 — (703O4 1 — qa4 1 — i/aia203a4 
Here we choose the order q < t < ai < 02 < • • • and we can not set i = 1 as this moment yet. 



(66) 
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Now it turns out to be expedient to start by eliminating 04. This can simply be done by omitting 
the factor 1/(1 — i/aia2a3a4) and making the substitution a4^i/aia2a3, obtaining 



G3iq,t) 



1 



1 — 9010203 1 — (70102 1 — (70103 1 — qai 
1 1 



1 



1 



1 — qt/ai 1 — qt/aias 1 — qt/aia2 1 — qt/aia2a3 

Setting i = 1 is valid here. Grouping terms containing the same subset of the variables 01,02, 03 gives 
1 1 

G3{q) 

i — qui 1 — y/ ui 

1 



1 — qai 1 — q/ai 
1 



1 — 90102 1 — 9/0102 

1 1 



(67) 



1 - 90103 1 - 9/0103 
1 



1 - 9010203 1 - 9/0102O3 
Likewise, we can easily see that the general form of ()66|) is 



Gkiq,t) 
with 



II l^qa,A{S)){^Jl,l-qA{, 



SC[2,fc] 



SQ[2,k] 



1 



qA{S)ak+i J I - t/aia2 ■ ■ ■ OfcO^+i 



A{S) = Y[a^. 



Removing the last factor and setting o^+i — </oi02 • • ■ o^ gives 



n i^qaiAis)) ( _n 1 



Gkiq,t) 

• SC[2,fc] 

and by setting t — I this can be rewritten as 



1 



SC[2,fe] 



qtA{S) / aia2 ■ ■ • Ofc 



Gkiq) 



lil-90iA(5) 1-9/01 



SC[2,A:] 



A{S) 



Now comes the next trick: grouping terms according as A{S) contains 02 or not. This gives 



Gkiq) 



n ^1 [ n - - 



aia2A{S) 



(68) 



a-,a"---a" 



To appreciate the significance of this step let us see what this gives for k — 3. Grouping terms in (j67|) 
as was done in ([55]) gives 



63(9) 



1 



1 



1 



1 



1 — 901 1 — 9/01 1 — 90103 1 — 9/01O3 

1 



1 



1 



1 



1 - 90102 1 - 9/0102 1 — 901O2O3 1 — 9/01O2O3 



(69) 



Let us now see what the partial fraction algorithm gives if we first eliminate 02. This entails 
computing the constant term 



Q = 



1 



1 



1 



1 



1 - 90102 1 — 9/0102 1 — 9010203 1 — 9/010203 
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Using the terminology of [3] we note that the first and third factors are contributing and the other 
two are duaUy contributing. Thus, 



Q 



A, 



1 - qaia2 1 - 5010203 



^1 +^3 



(70) 



ith 



A, 



9 9 



(0102 - q){l - gai0203)(oi0203 - q) 



1 



9 9 
0JO2O3 



(1 - gai02)(oi02 - 9)(oia203 - q) 
Using dTD]) in ^ gives 



a.=l/,a, {l-q'){l-a,){l-q^la,) 
_ 03 



1 



1 — qai 1 — q/ai 1 — (70103 1 — (7/0103 
111 1 



1 — (701 1 — (7/01 1 — (701O3 1 — (7/0103 



Ai + A3 
Ai + A3 



(71) 



The last equality is due to the fact that Ai and A3 do not contain oi. Next we will compute the 
constant term 

111 1 



1 — qai 1 — (7/01 1 — (70103 1 — (7/0103 



Q' = 



The surprise, which is the whole point of the factorization in (|68p . is that this leads to the same partial 
fraction decomposition! More precisely we see that 



Q' = 



Bi 



1 — qai 1 — (70103 



with 



Bi 



(01 - q){l - 50i03)(oi03 - q) 

^ ^ ^^1^3 

(1 - (70i)(oi - g)(oi0203 - g) 
Thus ([7T|) becomes 



= B1+B3 



1 



.,=1/, (1-(72)(1-03)(1-(7V«3) 

_ 03 

a. = lM,a3" («3-l)(l-'Z^«3)(l-g2) 



Al 
= A3. 





= A\ 


+ Al 


+ 2A1A3 










a% 





It is easy to see that the same collapse of terms occurs in the general case. Indeed we can rewrite 
in the form 



Gk{q) = 



n ^ ^ 



SC[3,fe] 



02^(5) 



J- J- l-qaiA{S) 1-g/oi 



SC[3,...,fc] 



A{S) 



We can see that, in both constant terms with respect to oi and 02, the first member of each pair of 
factors contributes and the second dually contributes, and the partial fraction algorithm yields 



n 

r" 1 — qaii 



1 



SC[3,...,fc] 



- 90102^(5) 1 — 9/0102^(5*) 



TC[3,...,fc] 
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with 



Ct = 1 - qaia2A{T) 



n 



1 



1 



SC[3,...,fc] 



n 



qaia2A{S) 1 - 17/0102^(5) 
1 1 



(1 - q/aia2A(T)) ^^J"^ 1 - 90103^(5) 1 - q/aia2A{S) 



a2 = l/qaiA{T) 



a2 = l/qaiA{T) 



(1 



n 

SC[3,...,fe] 



1 



1 



1 - A{S)/A{T) 1 - q^A{T)/A{S) 



and we see that, as in the case k = 3, all of these coefficients are independent of oi. Moreover we can 
also easily see that 



(1 - gaiA(T)) n 



1 



1 - qaiA{S) 1 - q/aiA{S) 



SC[3,...M] 

This reduces the computation of Gk{q) to the sum of 2*^^^ 



ai = l/qA{T) 



2 

j constant terms of the form 



2 E 

i<i<i<2'"-2 



AiAj 



Note that for fc = 5 we are reduced to the calculation of 2^ + (^2 ) ~ '^^ constant terms. Most 
importantly in each of these constant terms the denominators have at most 14 factors. The latest 
version of the partial fraction algorithm (whose update is motivated by the computation of G5{q)) 
posted in the web site 



h ttp:/ / www. combina tones, net.cn /homepage/xin / maple/el 12. rat 



computed these 36 constant terms on a Pentium 4 Windows system computer with a 3G Hz processor 
in about 22 minutes which is a considerable time reduction from the 2 hours and 15 minutes that took 
previous versions of the algorithm to compute these constant terms. 

The same approach can be used to calculate ^^5(9), but in a much simpler way. The constant terms 
have to be appropriately modified. Again we will start with the case fc = 3. 

The fc-tuple of sets interpretation of the constant term in (|26p given in Section [31 yields that to 
obtain Wk{q) we must compute the constant terms corresponding to the 2'' systems obtained by 
requiring each Ai to have 2 or more elements than its complement in all possible ways and then 
carry out an inclusion exclusion type alternating sum of the results. 

A moments reflection should reveal that to get W3{q) = W3{q, 1) we need only modify ([66)1 to 



Waiq, t) = ((1 - 04/01) (1 - 1/02) (1 - l/as)) 



1 



1 



1 



1 



1 — 9010203 1 — gai02 1 — 90103 1 — qai 

1 111 



1 



1 — 9020304 1 — 90204 1 — qa^an 1 — 904 1 — i/oia20304 
In fact expanding the first factor gives the 8 terms 



(72) 



a-(aAaAai 



1 - l/a2 - l/aa - a^/ai + 04/01(22 + 04/0103 + I/02O3 — 04/01O2O3. 
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And we see that the 8 constant terms obtained by expanding this factor in (|72p correspond in order 
to the foUowing 8 modified versions of S'^ 



l^ll 


= d 


l^il = d+1 


\Ai\ 


d 


\Ai\ 


d 


1^2 1 


= d 


\A2\ = d 


\A2\ 


= d+l 


1^2 1 


d 


1^3 1 


= d ' 


\A3\ = d ' 


1^3 1 


d ' 


1^3 1 


= d+l 




= d 


r^il = d-l 




d 




d 



\Ai\ 


= d + l 


\Ai\ 


= d+l 


\Ai\ = d 


l^ll 


= d4 


- 1 


\A2\ 


= d + l 


\A2\ 


d 


\A2\ - d+l 


1^2 1 


= d^ 


- 1 


l^3| 


d ' 


\A,\ 


= d + l ' 


IA3I - d + l ' 


l^3| 


= d^ 


- 1 




= d-l 




= d-l 






d- 


- 1 



Now the elimination of 04 in ([7^ and then setting t = 1 (as for 6*3(5)) gives 
H^3(g) = ((l-l/a?a2a3)(l-l/a2)(l-l/a3))x 



1 



1 



1 



1 



1 — (7010203 1 — goi02 1 — 90103 1 — goi 
1 1 1 



1 



1 - g/oi 1 - g/oi03 1 - 9/0102 1 - 9/010203 
For general fc, we are left to compute the constant term 

'11 



Wkiq) = (1 - l/o?02 ■ • • flfc) n (1 - 1/°^ 



n r 

SC[2,A:] 



90iyl(5) 1 -9/oiyl(5) 



Using this formula, the updated package will directly deliver ^5(9) in about 17 minutes. This is 
because the factors in the numerator nicely cancel some of the denominators of the intermediate 
rational functions. 
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